Solvable Stochastic Dealer Models for Financial Markets 
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We introduce solvable stochastic dealer models, which can reproduce basic empirical laws of 
financial markets such as the power law of price change. Starting from the simplest model that 
is almost equivalent to a Poisson random noise generator, the model becomes fairly realistic by 
adding only two effects, the self-modulation of transaction intervals and a forecasting tendency, 
which uses a moving average of the latest market price changes. Based on the present microscopic 
model of markets, we find a quantitative relation with market potential forces, which has recently 
been discovered in the study of market price modeling based on random walks. 



PACS numbers: 02.50.Ey Stochastic processes, 05.40.Jc Brownian motion, 
physics, financial markets, business and management 
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1. INTRODUCTION 

Research on financial markets using methods and con- 
cepts developed in physics has increased considerably 
over the last decade. Various kinds of stylized facts or 
empirical laws of markets have been discovered from high 
precision market data of gigantic size [H 3 01 Hi • The 
next goal of this econophysics study is to attempt to es- 
tablish the reasons for these empirical findings. Just as 
with the Boyle-Charles' macroscopic law which can be 
derived from a simple microscopic ideal-gas model, we 
hope to construct a simple microscopic model of a market 
that can reproduce major empirical findings. By relating 
macroscopic market behavior to microscopic dealers' ac- 
tions, we may find a pathway to control the markets, so 
as to avert bubbles and crashes, which occasionally cause 
problems in the market. 

The study of modeling dealers' action is carried out 
with so-called agent-based models. This approach is sup- 
ported not only by economists but also by information 
scientists and physicists Agent-based models 

can in practice reproduce dealers' actions in the mar- 
ket and they can also reflect empirical laws of markets to 
some extent. However, agent models generally include a 
huge number of parameters, and it has proved difficult 
to understand the relation between the parameters of the 
model and resulting market behavior. 

In order to find relationships between the parameters 
of dealers' actions and market behavior, we have already 
introduced a kind of minimal model of an agent-based 
market which consists of dealers with simple determinis- 
tic time evolution rules 0|[l3l[IH- With this model, we 
successfully reproduced most of the basic empirical laws 
using a minimal number of parameters, and found that 
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there are only three important effects needed to repro- 
duce the empirical laws. The first effect is the compro- 
mise pricing of both buyers and sellers, who tend to allow 
the particular transaction price they have in mind to ap- 
proach the current market price in order to make a deal. 
From this effect, transactions occur spontaneously in the 
market and the price rises and falls almost randomly. 
The second effect is the self-modulation of transaction 
intervals, that is, the rate of a dealer's clock depends on 
the latest moving average value of transaction intervals. 
When market activity becomes high, dealers accelerate 
their transaction rates, and by this effect we can repro- 
duce empirical statistical properties of transaction inter- 
vals which deviate from a simple Poisson process. The 
third is the trend- follow effect, that is, dealers forecast 
upcoming prices using the latest market trend which is 
defined by a moving average of price changes. This fore- 
casting effect makes the price change distribution follow 
a power law quite similar to that of the real market. 

In this paper we first introduce a stochastic version of 
the dealer model which is even simpler than the above 
(deterministic) model. In the case of the deterministic 
dealer model we needed at least three dealers to repro- 
duce market properties; however, in the present stochas- 
tic model we require only two. The advantages of this 
stochastic model are not only its simplicity but also its 
solvability by analytical calculation. In the usual agent- 
based approaches intensive numerical simulation is the 
only way to obtain results; in such cases exact or strict re- 
sults are rarely obtained. Based on this stochastic dealer 
model and its variants we can derive the major empirical 
results mentioned above, that have already been obtained 
by simulation of the deterministic dealer model by theo- 
retical analysis. 

Apart from agent-based modeling, the standard way 
to model markets is by utilizing random walks. It is now 
widely known that Bachelier introduced a random walk 
model for market prices five years earlier than Einstein's 
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random walk model of Brownian motion fisl . Work 
on portfolio theory, option price formulation IJ] and the 
ARCH and GARCH models [13 [H], which has led to No- 
bel prizes for their developers, are all based on random 
walk models. 

Recently, one of the authors (M.T.) has introduced a 
new type of extended random walk model of the market, 
the so-called PUCK model, in which a random walker 
moves according to a deforming potential force, the cen- 
ter point of which is given by the moving average of the 
random walker's traces. By using this generalization, all 
major empirical laws can be established; moreover, dy- 
namical behaviors such as bubbles, crashes and inflations 
can also be described as following from special cases of 
the market potential force l3| 3 ■ The ARCH model can 
also be derived as a special limiting case of this extended 
random walk modcl|19l|. 

Considering the wide applicability of the PUCK model 
has led to an open question concerning the origin of mar- 
ket potential forces. This question has been part ially 
answered by using the deterministic dealer model [20|. 
Here, we are able to provide quantitative answers by us- 
ing the stochastic dealer model. 

In the next section we introduce our stochastic dealer 
models step by step in sequential subsections. The third 
section is devoted to the relationship with the PUCK 
model, in which we will see how dealers' actions produce 
a market's potential force in a quantitative discussion. 
The final section contains a summary. 



2. THE STOCHASTIC DEALER MODEL 

In this section we introduce three stochastic dealer 
models, Modcl-1, Modcl-2 and Model-3. Model-1 is the 
simplest market model in which the framework of the 
stochastic dealer model is introduced. Then, we note 
two empirical properties which Model-1 cannot reflect. 
In Model-2 and Model-3 we introduce two additional ef- 
fects respectively to deal with these difficulties. After 
combining these revisions, the stochastic dealer model 
fully reflects all major empirical laws of markets. 



2.1. Model-1 

Firstly, we assume an artificial market consisting of 
only two dealers who are offering both buying and selling 
prices. The buying price, or bid price, is the current 
maximum price at which the dealer wants to buy. The 
selling price, or ask price, is the minimum price at which 
he will sell. For each dealer the ask price is always higher 
than the bid price, because they want some margin, and 
the difference between these prices is called the spread, 
which is assumed to be a constant, L, in this model. 
We define the i-th dealer's mid-price at time t, Pi{t), as 
the average of his bid and ask prices. When \pi{t) — 
P2{t)\ is less than L, these dealers do not transact as 



their transaction conditions are not fulfilled (FIGllJi). In 
such a case dealers are assumed to change their prices 
randomly according to the following rule. 




P^{t) + cMt) iei,2, 

Ap (prob. 1/2) 
Ap (prob. 1/2) 



(1) 



Here, fi{t) is a random noise for the i-th dealer, and c is 
a constant parameter. Then, the distance between pi{t) 
and P2it) is checked. If it is greater than or equal to L, 
then one dealer's bid price is higher than the other's ask 
price, and a transaction occurs (FIG[Tt). In such a case 
a unit volume deal is assumed to be made, with the mar- 
ket price given by the averaged price of the two dealers' 
mid-prices. After this transaction their mid-prices are 
assumed to shift to the market price. These processes 
arc repeated again and again and the time proceeds in 
unit of At. 

It should be noted that there is a possibility that this 
model produces a negative value for market price. In 
such a case the step width of the price change, Ap, should 
depend on the market price, such that the value of Ap 
is proportional to the market price to avoid crossing the 
origin. Here, we pay attention only to the case that the 
price fluctuation level is much smaller than the market 
price and for simplicity, we assume that Ap is constant. 





(c) (d) 
Transaction 



FIG. 1: Time evolution of the dealer model. Squares and cir- 
cles denote ask and bid prices, respectively. The i-th dealer's 
mid-price is denoted by Piit). (a): In this situation no trans- 
action occurs, (b): The dealers' prices follow random walks, 
(c): When the distance between pi{t) and P2{t) is greater 
than or equal to L, a transaction occurs and the market price 
is defined by the averaged price of the two mid-prices, (d): 
After this transaction both dealers' mid-prices move to the 
market price. These processes are repeated. 

For convenience of analysis we define another unit of 
time called the tick time, denoted by n, which takes an 
integer value incremented at each occurrence of a trans- 
action. Accordingly, P{n) denotes the market price at 
tick time n, and the n-th transaction interval, /(n), is 
defined by the time difference between the n — 1-th and 
n-th transactions. 
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In FIGlU we plot an example of resulting market prices 
and corresponding transaction intervals. In the sub- 
windows of these figures we also plot the probability den- 
sity function of price changes \P{n) — P{n—1)\ and trans- 
action intervals I{n) both on a semi-log scale. It is clear 
that the tail parts of both of these distributions are well 
characterized by exponential laws. As for the intervals, 
this result implies that the occurrence of transactions of 
this model is approximated by a Poisson process. It is 
interesting that the price change distribution of this sim- 
plest model follows an exponential distribution except 
aroimd AP = 0, instead of Gaussian distribution. 
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FIG. 2: Examples of time series of market prices (a) and 
transaction intervals (b). Sub- windows of these figures show 
the probability densities of market price changes and trans- 
action intervals on a semi-log scale. The parameters for this 
simulation are as follows: L = 0.01, c = 0.01, Ap = 0.01 and 
At = Ap- Ap. 

We can explain the functional form of the tails of these 
distributions as follows. We define the difference of the 
dealers' prices by D{t) = pi{t)—p2{t), then the condition 
for the occurrence of a transaction is given by > L, 

and we also define the mass center by the average of 
the two mid-prices, G{t) = {pi{t) + p2{t)} /2. As the 
mass center at the time of transaction gives the market 
price, market price statistics can be calculated from the 
information about G{t) for times > L. Namely, 

AP is given by AG which is defined by AG = G{t) ~ 
G{t') , where t' is the previous transaction time. These 
two variables define a two dimensional random walk with 
absorbing walls at D{t) = L and D{t) = —L as shown 
in FIGlHl The stochastic dynamics is described by the 
following set of equations 



D{t + At) 




(prob. 
(prob. 
(prob. 



1/4) 

1/2), 
1/4) 



(2a) 



AG{t + At) = AG(t) 




(prob. 1/4) 
(prob. 1/2). 
(prob. 1/4) 



(2b) 



When the random walker reaches one of the absorbing 
walls a transaction occurs, and by the transaction rule of 
Model-1 that the prices of the two dealers are then set to 
the market price, the random walker goes to the origin, 
and a new random walk begins. In this 2-dimensional 
formulation the transaction interval, /(n), is given by 
the survival time, that is, the time the random walker 
starting from the origin takes to reach one of the absorb- 
ing walls. Similarly, the market price change, AP(n), is 
given by the random walker's location on the G axis. 
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: AG 



FIG. 3: Random walk in the AG — D plane. The particle 
starts from the origin and continues a random walk until it 
touches one of the horizontal walls, meaning that a market 
transaction occurs. The transaction interval is given by the 
survival time of this random walker, and the market price 
change AP is given by the distance along the AG axis from 
the origin to the position of the particle on the absorbing wall. 

In Fig. 3 an example of a random walk is shown for 
better understanding of this mapping of Model-1 to a 2- 
dimensional random walk. Here, the horizontal axis is 
the X-axis and the vertical one is the y-axis. For theo- 
retical analysis we consider a continuum limit such that 
the mesh sizes of space and time go to zero. Under the 
condition At = (Ax)^ ~ (A?/)^, we find that the proba- 
bility density u{x^ y, t) of the particle in the [x, y) plane 
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at time t is described by the following diffusion equation 
du{x,y,t) 



u(x,2/,0) 
u{x, 0, t) - 



dt 

-- 6{x,y- L) 
u{x,2L,t) = 



4 dx^ dy^ 



:Initial condition 
:Boundary condition 



(3) 



(4) 



Here, is equivalent to the diffusion coefficient of this 
2-dimensional random walk. The initial condition is the 
delta function, and the boundary condition is given by 
the absorbing walls on the y axis, while there is no bound- 
ary in the x direction. This diffusion equation is solved 
exactly as follows, 



u{x,y,t) = 



1 



cL\/ nt 



n=l 



suY — s\i\Fny-e 



(5) 



61, •• 

m 



(3{k) is the Dirichlet beta function defined as 
(-1)" 

Using this result we can calculate 



E 



n=0 ' 



means and variances of transaction intervals and volatil- 
ities with results as shown in TABLE [D 

TABLE L Exact solutions for the means and variances of 
transaction intervals and the absolute value of price changes in 

1 9 

ModeH. K is Catalan's constant; K = 8(2) = - dO 

^ ' 2 Jo sine 





Average 


Variance 


Interval (/) 






Price change (|AP ) 


4:KL/-k'^ 





Here, Pn = We obtain the distributions of trans- 
action intervals Qi{I) and price changes (52(|AP|) by 
calculating distributions of survival times and absorbed 
points from ([5]). 



Qi{i) 

Q2{\AP\) 



^ II -1-1 



11=1 ^ 



. oo 

-y 



(-1) 



lAPI 



i^(6a) 
(6b) 



In the case of large values of / and |AP| in Eq.®, these 
summations are dominated by the term of n = 1. So, the 
functions / and |AP| can be approximated as 



Q\(l) oc e 



Q2(iAPi: 



= |AP| 



(7a) 
(7b) 



From these results we can derive the exponential laws 
of interval distributions and the price change distribu- 
tions already seen in FIGlH It is confirmed that the 
theoretical values of the decay constants (2L/c7r)^ and 
L/tt, fit well with the numerical results. 

Higher order moments of the distributions of trans- 
action intervals and price changes are also obtained ex- 
actly from Eq.®. The fc-th moments of < Q\(l) > and 
< Q§(|AP|) > are calculated as follows. 



< 



(/)> 



<g^(|AP|)> = 



1) 



r(2fc 
4L'=r(fc- 



E, 



rk+l 



(8a) 
(8b) 



Here ^{x) is the gamma function and Ek are 
the Euler numbers appearing in the expansion of 

secx = ^ ; £^0 = ^-.Ei = l,E2 = 5, £'3 

fc=0 



(2fc)! 



2.2. Model-2 

In this subsection we focus on the statistical differences 
between transaction intervals of real markets and those 
of Model-1. In real dollar- yen exchange market data pro- 
vided by EES for six years from 2000 to 2005, we find that 
the transaction intervals exhibit a circadian pattern even 
for Foreign Exchange markets which are open continu- 
ously as shown in FIG H] for the Dollar- Yen market. As 
can be clearly seen, large numbers of transactions occur 
during office hours of Tokyo, London and New York and 
the transaction density is least a little before the open of 
the Tokyo offices. 




10 12 14 

Time (hour) 



FIG. 4: Diurnal pattern of transactions in the Dollar- Yen 
market. The vertical axis depicts the mean number of trans- 
actions per quarter hour. 

In addition to this 24-hour pattern, there are fluctu- 
ations with much shorter time scales as typically shown 
in FIGini In this figure transaction events are shown by 
bars at the top and corresponding Dollar- Yen rates are 
plotted in the lower section. Here the window size is 
ten minutes and we can find places where bars tend to 
cluster, marked as "dense", and others where bars are 
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"sparse". The distribution of these intervals is clearly 
seen not to correspond to the simple theoretical model of 
a Poisson process. 
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FIG. 5: Dollar- Yen rates for ten minutes (lower curve), and 
transaction intervals (upper lines). There are dense and 
sparse periods for the occurrences of transactions. 

The clustering properties of transactions are known to 
be well modeled by a self-modulation process introduced 
by the authors [21|, which is described as follows. 

x{n + 1) = e{n)- < x{n) >r +f{n). (9) 

where r is the time scale of self-modulation, and 



N-l 



< x{n) >r= 



k), 



(10) 



A:=0 



e{n) and f{n) are independent noises, and N is the num- 
ber of transactions occurring within r seconds. This pro- 
cess is a modified Poisson process whose mean value is 
given by the moving average of transaction intervals over 
the past r seconds. As a result, there is a greater ten- 
dency to cluster and the so-called 1/f fluctuation is real- 
ized in general. 

Model-2 is designed to satisfy the real interval property 
by applying the self-modulation process. We estimate the 
distribution of e(n) from real data by using the following 
relation for the transaction intervals, I{n). 



e{n) 



I{n) 



< I{n) >r' 



(11) 



Here the typical value of r is 150 seconds. It is con- 
firmed from Dollar- Yen rate data that the distribution of 
e{n) follows an exponential distribution in general with 
mean value of unity. This exponential distribution is 
favorable for our model construction as Model- 1 auto- 
matically produces the exponential interval distribution. 
As we can control the speed of transaction intervals by 
controlling the speed of diffusion, we obtain a revised 
model, Model-2, by modifying the constant parameter c 
in Eq.ll]), which is directly related to the diffusion co- 
efficient, making it a time-dependent parameter c(n) as 
follows. 




p,{t)+c{n)f,{t) iel,2, 

+Ap (prob. 1/2) 
-Ap (prob. 1/2) 



(12) 



where 



c(n) 



< / >c=l 
< I >r 



(13) 



In Eq. (|13p . < / >c=i is a mean of transaction in- 
tervals shown in TABLE [T] in the case of c = 1 and 

< / >c=i= i^/2. < / >^ is the moving average of 
transaction intervals averaged over the latest r seconds 

defined as < / >r= — -^(^t- ~ k). In this equation, 

I{n — k) is the transaction interval that is the fc-th tick 
earlier than the n-th tick. TV is the number of transac- 
tions within T seconds from time n. If I{n) > t, we set 

< / >r= I{n)- It is known from Eq. p^ that for larger 

< / >c=i over < / >T, the value of c{n) is larger, that 
is, dealers tend to make larger changes their prices to ef- 
fect more rapid transactions when transaction intervals 
become shorter in the market. By this effect, Eq.([2]) of 
Model- 1 is modified to 



D(t + At) = D(t) 



AG{t + At) = AG{t) 



~2c{n)Ap (prob. 1/4) 
±0 (prob. 1/2), (14a) 

-2c(n)Ap (prob. 1/4) 

~c{n)Ap (prob. 1/4) 
±0 (prob. 1/2). (14b) 

-c{n)Ap (prob. 1/4) 



Examples of random walk traces are shown in FIG (HI 
As known from this figure the initial condition and the 
boundary conditions are the same; however, the step 
size changes for each random walk following the self- 
modulation formulation. By this effect the transaction 
intervals tend to form clusters as shown in FIGl?] 

In order to make the interval distribution fit well with 
that of the real Dollar- Yen market we introduce two 
thresholds for the value of < / When < / >r< 3, 
we set < / >T= 3, and when < / >t> 50, we set 
< / >r= 50. These restrictions are needed to prevent 
intervals from converging to zero, or from diverging to 
infinity. With these minor revisions an example time 
sequence produced by Model-2 is plotted together with 
one produced by Model-1 in FIGlT] Comparing these 
two sequences, we observe that Model-2 can reproduce 
the clustering property quite well. Moreover, the lower 
interval sequence looks similar to the real sequence shown 
in FigO Actually, the distribution of transaction inter- 
vals arising from Model-2 is now very close to that of the 
actual interval distribution as shown in FIG151 
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FIG. 6: Examples of random walk traces arising from Model- 
2. Compared with FIG |3] for Model- 1, the step size depends 
on past transaction intervals. 
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FIG. 8: Comparison of the distributions of transaction inter- 
vals, (a) An example of a time series of intervals produced 
by Model-2. (b) Cumulative distributions of transaction in- 
tervals on a semi-log scale. Actual data plotted by dashed 
line during 9:00-10:00, New York time. The distribution of 
Model-2 plotted by the solid line has a longer tail than the 
exponential distribution produced by Model-1 in the dotted 
line. The distribution of real intervals represented by the 
dashed line is close to that of Model-2 with the parameters 
L = 0.01, r = 150, Ap = 0.01 and At = Ap ■ Ap. 



FIG. 7: Examples of transaction occurrences. The upper se- 
quence is produced by a simulation of Model-1 and the lower 
one by Model-2. 



2.3. Model-3 



In this subsection we shift our attention from transac- 
tion intervals to price changes. We know that the price 
change distribution of Model-1 is characterized by an ex- 
ponential distribution while that of the real market is of- 
ten characterized by a power law. It has been established 
that such power law distributions can be derived by intro- 
ducing the effect of trend- following prediction This 
effect can be introduced to our stochastic dealer model 
by simply adding a further term, d < AP >j\/ At, which 
is defined as follows. 



P^{t + At) 

m 

where 



+ d<AP >M At + cf,{t), (15) 
j+Ap (prob. 1/2) 



-Ap (prob. 1/2) 



«e 1,2, 



M~l 



— (^^ - k)AP{n - k). (16) 



fe=0 



Here AP{n) = P{n) — P{n — 1) is the price change at the 
71-th tick. The new term, < AP >a/, is a kind of moving 
average of price changes for M ticks with weights that 
decay linearly. The parameter d in (jlSp is an important 
parameter that governs the dealers' strategy. A dealer 
with positive d is a trend-follower who predicts upcoming 
market prices proportional to the latest price slope. On 
the other hand, a dealer with a negative d is called a 
contrarian who forecasts that upcoming market prices 
will go against the trend and that the present market 
price is close to a local maximum or minimum. 
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Adding this effect, the equations ([2]) in Model- 1 are 
modified as 



D{t + At) = D{t) 




(prob. 1/4) 
(prob. 1/2) , (17a) 
(prob. 1/4) 



AG'(i + A<) = AG(i) + d < AP >m At 
-cAp (prob. 1/4) 
±0 (prob. 1/2) . (17b) 
-cAp (prob. 1/4) 

In the 2-dimensional random walk representation the ini- 
tial conditions and the boundary conditions are invari- 
ant, however, we have a horizontal flow proportional to 
d < AP > M as shown in FiglQl The existence of this flow 
implies that the distance of the absorption point from the 
origin is greater than for the original Model-1. As the 
vertical motions are completely identical the transaction 
interval is also identical; however, the absorbed point on 
the horizontal axis is shifted by I{n)d < AP >m- The 
strength of the flow depends on the parameter d and the 
latest price changes. 




AG 



FIG. 9: Random walks on AG - D space for Model-3. Two 
traces are plotted for comparison with the results of Model-1: 
The dotted line is for Model-1 and the solid line is for Model-3, 
both produced using the same random number generator. 

In this revised model the transaction intervals are iden- 
tical to those of Model-1 because Eq. (|17ap is the same as 
Ea. (Pa|) . while the market price change is described by 
the following equation. 



AP{n + 1) = I{n) ■d<AP>M +F{n). 



(18) 



Here, the first term of the right hand side is the distance 
covered by the flow, and d < AP >m gives the intensity 



of the flow, I{n) is the transaction interval. The second 
term is identical to the price change of Model-1. From 
the results already obtained for Model-1 it is clear that 
both I{n) and F(n) are random variables characterized 
by exponential functions, so Eq. (|18[) follows a random 
multiplicative process. We know that a time series which 
is produced by a random multiplicative process generally 
follows a power law if the process satisfies a stationary 
condition. In particular, in the case that M = 1 and 
I{n) and F(n) are independent in Eq. p^ . we have an 
exact solutionjl^l. In this Model-3, I{n) and F{n) are 
not statistically independent and so we do not have an 
exact solution; however, if the effect of F is negligibly 
small, then we can approximate the result using the so- 
lution for the independent case. We have the exponent /? 
for the power law of cumulative price change distribution 
as follows. 



\d\f^ <I{n)^ >= 1. 



(19) 



Here < I{n)^ > is the /3-th order moment of I{n), so we 
can apply Ea. ([8a|) to Eg. pQ)) . As a result, we have 



' ' c2/3r(2/3 + l)" 



= 1. 



(20) 



It is known that empirical values of the power exponent 
of the cumulative price change distribution are around 
—3 in the actual market, so we set /3 = 3,L = 0.01 and 
c = 0.01 in Eq.jlQl). Then we have \d\ - 1.25. We can 
reproduce the power law of price change distribution with 
exponent —3 as represented in FIGfTUl Other parameters 
are taken to be Ap = 0.01, At = Ap ■ Ap and M = 1. 



3. RELATION BETWEEN THE DEALER 
MODEL AND THE PUCK MODEL 

We have seen that the stochastic dealer models can re- 
produce important empirical features of markets. In this 
section we examine the relation to the market potential 
T3| 3- In our previous work, we 



model called PUCK 
showed that the deterministic dealer model can repro- 
duce market potentials confirmed by numerical simula- 
tion, and the essence for reconstruction of the potential 
is found to be the dealers' forecasting effect using moving 
averages [l^. We now present an analysis based on the 
present stochastic dealer model. 

It is easy to confirm that price changes produced by 
Model-3 yield non-trivial market potential functions as 
shown in FIG IllI In this figure the parameter d in 
Eg. p^ is changed. We set the parameter d = —1.0 
during the period from n=l to 1000 ticks, that is, the 
dealers' are contrarians who predict that the future price 
will move against the latest trend. The parameter d = 
in the period n=1001 to 2000 ticks, that is, the dealers 
arc simple random walkers. And d = 1.0 during 2001 
to 3000 ticks, that is, the dealers are trend-followers pre- 
dicting that near future prices will be proportional to 
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FIG. 10: Price changes produced by Model-3. (a) An exam- 
ple of time evolution of market prices produced by Model-3. 
(b) Cumulative distributions of price changes using a log-log 
scale. The dotted line is for Model-1; the solid line is for 
Model-3 and the dashed line for real data. The slope -3 is 
shown as a guide-line. 



a moving average of price changes. We can clearly ob- 
serve a stable, a flat and an unstable potential function, 
respectively, as expected. During 3001 to 4000, we set 
d = 1.0 when the average of the past M price changes, 
< AP >M, is greater than or equal to zero, and we set 
d = —1.0 when < AP >m< 0. In such an asymmetric 
case, we can find an asymmetric potential as shown in 
FIGfTTH: in that case the market price increases nearly 
linearly on a large scale. It is apparent that the market 
potential function and the dealers' forecasting effect are 
also deeply related in this stochastic model. 
The PUCK model is formulated as follows: 



P{n+1) = P(n) 



d 

■^U{P) \p=P{n)-PM + l{n) 



U{P) 



b{n) 
2i\7 



p2 



-F'{n), 
(21) 

(22) 



Here P(?t.) is the noise reduction price introduced by 
Ohnisi et. aL, also referred to as the optimal moving 
average price[23|, U{P) is the potential function defined 
by Eq. (P^ . F'{n) is an uncorrelated noise term, and 
PM+i{n) is the simple moving average over M + 1 ticks: 



PM+i{n) 



M 



1 ^ 

+ 1 ^ ^ 

fc=0 



k). If the market poten- 



tial is asymmetric as in FIG lTlJ d). we define the poten- 
tial function for a; < and x > respectively by using 
quadratic functions. In the case of symmetric potential 
as FIGfTTbbc. We can transform ([2T|) to the following 
equation: 



AP(n) 



bin) 



< AP >M +F'{n). 



(23) 



< AP >M is defined by Eq.([Tni). We note that Eq.^H]) 
in Model-3 and Ea. (P5)) have the same form of a linear 
stochastic equation, so the statistical property is indepen- 
dent of the noise property. Comparing the coefficients of 

< AP >]\f in Eq.llHI) and Eq.([23]), we have the simple 
relation: 

b{n) = -2d- I{n). (24) 
By taking the average over tick times n, we have 



< 6 >= -2d- < I >= -d\ - 



(25) 



where < a; > denotes the average of x over tick time n. 
This equation implies that the market is unstable (b < 0) 
when dealers are trend followers {d > 0) while the market 
is stable when dealers are contrarians (d < 0) . This result 
is consistent with our previous simulation results using 
the deterministic dealer model. 

In the PUCK model, when the potential is of quadratic 
type such as in FIGIllb.bc. the diffusion coefficient of 
market prices, a, is theoretically calculated as a func- 
tion of the potential coefficient 6[2J]. By introducing the 
newly derived relation. Eq. ((26|) . into the formula, we have 
a theoretical evaluation of the price diffusion coefficient 
for our Model-3. 



cr<j(A7l) 



2c2 



2c2 - dL2 



crd=o(An). 



(26) 



From this relation we find that in the range of \d\ < 
2c^/L'^ the value of a is finite, and we can expect the 
market price to follow a normal random walk over a long 
time scale. When d > 2c^/L^ the above formula is mean- 
ingless; however, the actual market price moves nearly 
monotonically or even exponentially as occurs for bubbles 
in real markets. Actually we can generate a bubble- like 
phenomenon by setting d > 2c^/L^ as shown in FIGfT2l 
Here, the time evolution is well approximated by an ex- 
ponential function as predicted by the PUCK model [2^. 
We do not investigate this phenomenon further in this 
paper, but it should be stressed that our model can de- 
scribe not only normal states of markets, but also ab- 
normal states such as bubbles and crashes where prices 
move monotonically. by tuning our model's parameters. 
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4. DISCUSSION 

We have introduced a new stochastic dealer model 
which consists of only two dealers, and showed that ba- 
sic empirical laws of financial markets are well reflected 
in terms of transaction intervals and price changes. In 
Model-1, both dealers change their prices randomly, so 
fluctuations of transaction intervals and price changes are 
also random. We calculated these statistical properties 
exactly. The occurrence of transactions is well approxi- 
mated by a Poisson process, and the price change distri- 
bution is well described by an exponential distribution. 

In order to make our model more realistic the follow- 
ing two feedback effects were introduced. One was the 
feedback effect of transaction intervals, and the other was 
the feedback effect of price changes, both caused by the 
dealers' observations of the latest market status. As a 
result, the basic model's random noise properties were 
modulated and the artificial market reproduced both the 
distributions of transaction intervals and price changes 
to follow long-tailed distributions which are quite similar 
to those of real markets. 

Moreover, we established that the dealers' action of 
prediction by using a moving average of past price 
changes generates market potentials in the PUCK formu- 
lation, and we derived a simple theoretical relation be- 



tween the stochastic dealer model and the PUCK model. 
Namely, we found the relation between the microscopic 
dealers' strategy and the macroscopic market's stability 
as defined by the PUCK model. 

As an application of the relation to the PUCK 
model, we checked the condition when market prices in 
our stochastic dealer model exhibit a bubble-like phe- 
nomenon in which price motion is approximated by an 
exponential function rather than a random walk. This 
transition from a random walk phase to an exponential 
growth phase is considered to be quite useful for discus- 
sions concerning how to realize a stable market. 

We expect that our models can be used as a base for 
market experiments. For example, we may be able to 
observe the effect of governmental intervention by intro- 
ducing a third dealer who only buys dollars over a given 
period. From the viewpoint of numerical simulation it 
is quite easy to increase the number of dealers, each of 
whom will have his own strategy. So we can construct any 
experimental market by adding or subtracting specific 
dealers, and observe the change in the market's macro- 
scopic behavior as a result. For example, we may be able 
to find a way to avoid market crashes by introducing a 
specially designed dealer who acts to stabilize the mar- 
ket. This kind of market experiment may contribute to 
future attempts at real market stabilization. 
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FIG. 11: Market potentials estimated for the time series produced by our stochastic dealer Model-3. The parameter values are 
changed: d = -1.0 from 1 to 1001 ticks, rf = from 1001 to 2000 ticks, d = 1.0 from 2001 to 3000 ticks, and from 3001 to 4000 
ticks we set d = 1.0 if < AP >m> and d = —1.0 if < AP >a/< 0. Other parameters are given as follows: c = 0.01, L — 0.01, 
Ap = 0.01, At — Ap ■ Ap and M = 10. To estimate the market potential function, we use 500 ticks and Af = 10 and we show 
market potentials for the periods (a), (b), (c) and (d). The dashed line is approximated by a quadratic function, y = ax^. Over 
the period (d), we fit the line in a; > and a: < respectively by using quadratic functions as well as in (a), (b), (c). 
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FIG. 12: A bubble-like phenomenon produced by Model- 
3. The dashed line is an exponential function of time as 
a guideline, y = exp(0.004a;) + 99. The parameters satisfy 
d > 2c^/L^; d = 2.0, c = 0.01 and L = 0.01, Ap = 0.01, 
At = Ap- Ap and M = 10. 



